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A classical particle oscillating in an arbitrary high-frequency or static field effectively exhibits a 
modified rest mass nicft derived from the particle averaged Lagrangian. Relativistic ponderomotive 
and diamagnetic forces, as well as magnetic drifts, are obtained from the niaH dependence on the 
guiding center location and velocity. The effective mass is not necessarily positive and can result in 
backward acceleration when an additional perturbation force is applied. As an example, adiabatic 
dynamics with mn > and mn < is demonstrated for a wave-driven particle along a dc magnetic 
field, my being the effective longitudinal mass derived from niaft. Multiple energy states are realized 
in this case, yielding up to three branches of my for a given magnetic moment and parallel velocity. 
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I. INTRODUCTION 

A large number of problems connected with multi- 
scale adiabatic dynamics of classical particles in oscilla- 
tory and static fields enjoy critical simplification within 
the guiding-center approach, which allows separating fast 
oscillatory motion of the particles from their slow trans- 
lational motion [T1|21|31I1|31|51|71[H1[S]- Often, the 
average forces on a guiding center are then written in 
terms of effective potentials such as ponderomotive 
[Ulinilllllllj, diamagnetic [TF, or others P 15]. Yet the 
applicability of the potential approximation is limited to, 
at most, nonrelativistic interactions, another drawback 
being the unphysical difference in fictitious fields — V^* 
seen by different species. 

The purpose of this paper is to offer an alternative 
approach that allows embedding the average forces into 
the guiding center properties, through redefining the par- 
ticle rest mass moff. The average acceleration is then 
attributed to the effective mass variations, which are 
naturally different for different species; hence no ficti- 
tious fields are introduced. By definition, this "object- 
oriented" formulation [14) is also intrinsically relativis- 
tic. Therefore, it equally holds for arbitrary adiabatic 
interactions [15], thus proving to be more fundamental 
as compared to the effective-potential approach. 

Previously, the effective mass TOeff was similarly in- 
troduced for an electron interacting with an intense laser 
wave in vacuum, with additional fields considered only as 
perturbations [Tl[T71[Ta[Tl[2SllIllll[131[Mll2Sl[2i[271 
UHllin]- In this paper, we show that meff can be defined 
as well for any other multiscale dynamics of a particle 
in high-frequency or static fields. We offer a general for- 
mula for the effective mass and show how manipulations 
of THeS as a function of the guiding center variables yield 
the average forces and particle trajectories. We also show 
that the effective mass is not necessarily positive and can 
result in backward acceleration when an additional force 
is applied. As an example, we explore the average motion 
of a laser-driven particle immersed in a dc magnetic field. 
Multiple energy states are realized in this case and yield 
up to three branches of rrics and the effective longitudi- 



nal mass my for a given magnetic moment and parallel 
velocity. We show that both mn > and my < are 
possible then, the negative-mass regime too allowing for 
adiabatic dynamics. 

The paper is organized as follows. In Sec. |TT] we de- 
rive the general formula for mcff and the guiding cen- 
ter Hamiltonian accounting for additional perturbation 
fields, if any. In Sec. |III| we apply the effective mass 
formalism to the particle motion in a static magnetic 
field and rederive both the particle Hamiltonian and the 
magnetic drifts. In Sec. IV we explore the average mo- 



tion of a laser-driven particle in a static magnetic field 
and demonstrate the possibility for adiabatic dynamics 
at negative rUeS and my. In Sec. [Vj we summarize our 
main ideas. Supplementary calculations are given in Ap- 
pendixes: In Appendix |A] we obtain the general form of 
the drift Lagrangian employed in Sec.|TT] In Appendix [B| 
we show how the effective mass formalism allows deriva- 
tion of ponderomotive forces in various cases of interest. 



II. EFFECTIVE MASS 

Consider a particle undergoing arbitrary quasi-periodic 
oscillations superimposed on the average motion. In the 
adiabatic regime, one can map out the quiver dynamics 
by changing variables [I1[1I1I1[31IS1IZ1[H]; hence the 
guiding center is treated as a "dressed" , or quasi-particle. 
Suppose, for now, that the background fields causing the 
oscillations do not vary along the trajectory. Then the 
associated field tensor F^i, will not enter the averaged 
equations as a force. However, it will affect the motion 
such that, in response to additional perturbation fields 
Ffii,, the guiding center will react as if it had a modified 
mass. 

The effect is shown as follows. At zero F^i,, the average 
dynamics is determined only by the field tensor i^^^^ seen 
by the particle in the guiding-center rest frame K' , fur- 
ther denoted by prime. The action increment dS — Cdt 
in the laboratory frame K is then written as 



dS{A^,w)=dS'{F'^^) + dg, 



(1) 
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where is the four-potential such that Ff^i, = d^A^ — 
di^Afj^, V = (v) is the guiding center velocity in K, v is the 
particle true velocity, (...) denotes the time-average, and 
dG depends on the selected gauge. Use dt = ^ dt', where 
dt' is the time interval in K', 7 = (1 — u^/c^)"^/^, and c 
is the speed of light. Then the guiding center Lagrangian 
reads 



(2) 



where £' = dS' /dt' . Omitting an insignificant full time 
derivative, one can rewrite Eq. Q as [SOj 



(3) 



Hence L is formally equivalent to the Lagrangian of a 
free particle with an effective mass 



(4) 



that is, by definition, both gauge- and Lorentz-invariant. 
By definition, CJ — C'{F'^^), where F'^^ can be written 



in terms of A'^; thus 



Ca,, - Ga, 



(5) 



v/1 -«Vc2' 

for any v. Consider t) 0; in this case, A^ — 

C'iF'^,)^[CA'^^-{dGA'Jdt)]^^^, (6) 



so 



where we removed the subindex "A^" in the left-hand 
side, as C is gauge-invariant. 

In the absence of oscillations, C must equal — mc^, 
where m is the true mass; therefore, 

C = -mc'^ + CA'^{v^O)-CA'^{v^O). (7) 

For clarity, we assume below that L{v = 0) = ~mc^ . 
Then, using Eq. (A6), one can write TOoff as 



TOcff = ^ ( J • 1/ - (L) 



(8) 



where the right-hand side is to be evaluated in K' (hence 
the index " w = 0" ) ; J are the actions and u are the fre- 
quencies of oscillations in canonical angles, if any, to aver- 
age over (Appendix [A| . Therefore, apart from the latter, 
rrieff is proportional to the gauge-independent part of the 
averaged Lagrangian in the guiding- center rest frame |31j . 
Since C is calculated in i(r'(v), TOos is generally a func- 
tion of the velocity v. When slowly varies with the 
guiding center coordinate r or time t, rrics may similarly 
depend on those as well, so Eq. ([3| will automatically 
yield the average forces [Eq. ( |A8| ]. 

Suppose now that a particle interacts with a perturba- 
tion field F^i/ governed by = (A, (p), which is imposed 
over F^i, (TBI 1231 [12] • In the adiabatic regime, the orbit 
is not altered on the oscillation time scale; thus. 



-niesc^ ^/l - u2/c2 -t- - (v • A) 



(9) 



(e being the particle charge), and a nonelectromagnetic 
potential can be added similarly. Then, the canonical 
momentum equals P = p + (e/c)A, and the kinetic mo- 
mentum p is given by 



p = 7meffv 



(10) 



7 9v 

Correspondingly, the Hamiltonian H = P ■ v — C reads 



U = 7meffC - — V • -—- 
7 V ov 



eip, 



(11) 



and £ = 7i(f,P,i) is conserved when Ti is independent 
of time. Thus, meff can be viewed also as the normalized 
quasi-energy of an unperturbed [F^^, = 0) particle in the 
guiding-center rest frame, TOoff — £'!<?■ 



III. STATIC MAGNETIC FIELD 

Let us demonstrate how the effective mass formalism 
applies to the problem of particle motion in a dc magnetic 
field B = V X A, where 



TOc' e , ^ , 
L= + - v-A , 

7 c 



(12) 



and 7 = (1 — I . Assuming a smooth B, the 
motion can be averaged over Larmor oscillations at fre- 
quency f2 = eB/mc^, so the guiding center dynamics 
that remains is one-dimensional, and the associated ac- 
tion J = (mc/e) fi is conserved |13j . where jj, = p]_/2Bm 
is the magnetic moment, and p_L = 7mv_|_ is the rela- 
tivistic kinetic momentum transverse to B. Thus, the 
effective mass reads rrics = [fiB/j' — {L)']/c^, where the 
prime denotes the guiding-center rest frame K' ; /i and B 
are Lorentz invariants, and 



7' = Vl + 2/iB/mc2 (13) 

is constant. Since {L)' — —{mc^ + fiB)/^, one obtains 

nics = m^/l + 2fiB/mc'^, (14) 

which is a relativistic invariant. The guiding center mo- 
mentum is then given by jmcffv — Pp where we used 
a Lorentz transformation 7 — 77'. Hence the Hamilto- 
nian (fTT| reads |Ml [33] 



n = ^m^d^ + 2fj,Bmc'^ + p^c^ + e(p, (15) 

yielding, after omitting an insignificant constant, the well 
known nonrelativistic limit 



n ^ ^ p\ + [iB + e(p. 
2m " 



(16) 



A more precise calculation also delivers particle drifts 
pn, 1131(51] : Allow arbitrary Wy, yet assume nonrelativistic 
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vj_ so as to treat the transverse drift as a perturbation. 
Following Ref. |35j . we write the new guiding center La- 
grangian as £ = Cq + C\nt, where Lq — — moffC^/7 is that 
for the unperturbed motion, and 



Ant = -(v • A) - 



(17) 



is the interaction Lagrangian small compared to Cq. For 
simplicity, assume static fields, so the guiding center 
quasi-energy (111 is conserved. In this case, we can con- 
sider V as a function of f ; hence 5S = Q yields 

2 

(v -VJp = - V X B - ^ Vmcff, (18) 
c 7 

where V differentiates with respect to r at fixed v, and 
V* = V + Y^ii^ Vi){d / dvi) is the full spatial derivative. 



Eq. ( 18 ) is equivalent to 

vxB* = 0, B* = B + (c/e) V*x p, 
which can as well be put in the form 

_ _ B + (c/e)V,x(p||b) 
B + (c/e)p|, b-(Vx b) 



(19) 



(20) 



where b = B/_B, and p k, p^h. This generalizes a sim- 
ilar analysis, which was proposed in Ref. [35 for w <C c, 
to any v, such that v±_ <^ c. 

Eqs. ( 19 1, (20 1 yield the known expressions of the tra- 



ditional drift approximation [1 [331 [Ml ISZl [Ml l3Sj • How- 
ever, they also allow for an arbitrary rrieff, not necessar- 



ily that given by Eq. (14 1; thus additional strong fields, 
if any, are as well embedded here. Derivation of time- 
dependent and fully relativistic magnetic drifts [101 [S] 
using the effective mass formalism should be possible, 
too, but remains out of the scope of the present paper. 



IV. RELATIVISTIC WAVE FIELD OVER A 
STATIC MAGNETIC FIELD 

The unified effective mass formulation readily yields 
the ponderomotive forces previously derived from other 
considerations (Appendix [b| . In this section, we contem- 
plate another example of particle ponderomotive dynam- 
ics, which exhibits unusual properties that, to our knowl- 
edge, have not been covered in literature. 



A. Basic equations 

Consider a relativistic particle in a wave propagating 
along a static magnetic field j42l|43l|44l[45l|4Sl[471[48l[49l 
I5D] . Assume a smooth magnetic field B = V x Adc, ap- 
proximately in the z direction; then the vector potential 
Adc can be considered a linear function of the particle 
displacement from the guiding center location: 



1 



Adc = 7.B{z) (z X r_L). 



(21) 



For simplicity, we will also assume a vacuum electromag- 
netic wave with circular polarization in the plane trans- 
verse to B, so the corresponding vector potential reads 



mc 

e 



Qq 
V2 



(x cos ^ - y sin , 



(22) 



where the invariant Cq = eEQ/mcLu is allowed to slowly 
vary in space and time, Eq is the amplitude of the electric 
field E — — (l/c)(9Aw/<9i), and — cut — kz is the phase, 
with k — Lu/c. In this case, the particle motion is fully 
integrable, and the problem can be solved analytically. 

According to Sec. [IT] we calculate rrieff for uniform 
fields, using Eq. (12 1 with A = Adc + Aw. The effective 



mass is determined by the averaged Lagrangian in the 
guiding-center rest frame K' (further denoted by prime), 
which is found as follows. Since A depends on z and t 
only via ^{z,t), there exists an integral 

u^^~Pz/mc, (23) 

yielding that the following equality holds for any / [TB] : 

(/)-(7/)c/(7)c. (24) 

(The subindex ^ denotes that the averaging is performed 
over the phase rather than time.) Take / = L'; then 



1 



m(? — 



mc 



(P' • A'), 



With / V, Eq. (241 also yields 



(p)./m(7),, (p'), = 0, (7% 



(25) 



(26) 



Given that the average motion is solely in the z direction, 
K' is now defined as the frame where (p^)^ = 0. 

Hence the particle motion can be written as follows: 



xq + TZ cos 6 



sin^ 



y = Uo -TZsinO 



kuV2 (1-cr)' 
ao cos ^ 



zo 



ku 



Px = —V sin 6 — 



kuV2 (1-ct)' 
TZao o'cos(^ — 9) 

(1-^)' ' 

mcao cos ^ 



Py 



Pz 



-V cos I 



mcpQ 



V2 (l-aY 
mcttQ sin ^ 

Vao sin(C - 9) 
u^/2 (1-cr) 



(27a) 
(27b) 
(27c) 
(27d) 
(27e) 
(27f) 



Here po — uv/{c — v) is the normalized phase-averaged 
momentum; TZ = V/mflo, V, and 9 = a^ + 9o denote the 
gyroradius, the transverse momentum, and the phase of 
the free gyromotion superimposed on the wave-induced 
oscillations, a being a Lorentz invariant: 



no 

LOU 



^^o/7 



UJ - 



kv. 



(28) 
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FIG. 1: (color online) Effective mass rrieff of a wave-driven 
particle in a magnetic field [Eq. (|31[); m units] vs. oo and a. 



fi = 7'^/2_Bto is the invariant action of this gyromotion 
conserved under adiabatic perturbations [cf. Eq. (B6l]; 
Xq, Ho, zq, V, V, Oq are determined by initial conditions. 

To find u, substitute ^ — u + pz into 7^ = 1 + + p^, 
where p = p/mc, and average over ^ using (pz)^ = po- 



{u + pof = 1 + pI + {pI) 



Then, 



h = 



(29) 



(30) 



2(l-a)2' 

where a — cr{u), = 2p,B/mc^ is an invariant, and 
^mc^aQ equals the zero-B nonrelativistic ponderomotive 
potential $ = e'^ / Amuj'^ Eq. (|B3|]. 



Combining Eqs. p5^-(30l, one gets 



"T-cff 



a2(2- 



4(1- (7) 



2(1 -a)2 



(31) 



which is a covariant form of mes for the effective mass 



is expressed as a function of Lorentz invariants. Eq. (31 ) 



yields Eqs. (14 l-(|20 | at oq = 0, Eqs. (|B13|, (|B14f at B = 
0, Eqs. (|B9|, |b10P at i;/c < 1, and Eqs^ p3f (|B5| at 

0. 



v/c <ti I together with B — 0. From Eq. (|31|, it also read- 
ily follows that mcff(ao < 4) > 0; yet mcff(ao > 4) < at 
least for some a > 1 (Fig. [T|). Other properties of rrics 
are discussed in Sec. IIV Bl and C. 



B. Tristability 

In the presence of both relativistic effects and nonzero 
B, the cyclotron resonance is essentially nonlinear and 
permits multiple energy states at given v and p. To see 
this, rewrite Eq. (30 1 as 



U{u) EE 2{u-aofih 



1) 



2 2 



0, (32) 



s 3 



3 



FIG. 2: 



(color online) Solution of Eq. ( 32 1 for u vs. ctq = 
at (To 00, one has u\^2 ~ cro ± aoh/y/2, and 



00, one has ui,2 ~ cro ± 
~ /iVl + s^, where Ui correspond to different branches; 
it(0) = h^l + -I- ao/2. For a given ao, «2 and 113 appear 
simultaneously behind the nonlinear reso nan ce u — ao (gray, 
dashed) ; the condition is yielded by Eq. ( 33 1 . Shown here is 
the case ii = (so u = {7)^), s = 0: ao 
(green), and ao = 3 (blue). 



0.3 (red), ao = 1.5 



where ctq = ilo/^ |i51j. Eq. (32) is a fourth-order alge- 
braic equation; thus it allows up to four values of u, which 
also can be found analytically [51] §1.8-5]. (Explicit so- 
lutions are not shown here because of their complexity.) 
Since U{0) < and W(±oo) = +00, two solutions always 
exist, one of them being unphysical {u^ < 0). Further 
consideration of the signs in Eq. (32 1 shows [5H §1.6-6(c)] 



that, apart from degenerate cases, there exist either one 
or three positive roots, Ui > U2 > u^- Therefore, one or 
three energy states are possible (Fig. [2| , allowing for hys- 
teretic effects [53 IMl , which also have quantum 

analogies in solid-state physics [ST j l58 l l59] . 
The condition for multiple Ui reads 



(To > CTc = /i 



(l + s2)l/3 + (a2/2)V3 



3/2 



(33) 



Therefore, assuming that three branches exist for a given 
s, there must exist the same number of energy states for 
s = 0. The latter energy states correspond to three equi- 
libria in the momentum space {p±cos'>p, p± sin ip, pz), 
where ip — ^ + X being the gyrophase. Apart from 
the degenerate case when U2 approaches U3, the particle 
trajectory (27) is a continuous function of the initial con- 
ditions for each branch. Thus, assuming negligible dissi- 



pation (Sec. IV C I, all the three equilibria are stable here 
(Fig. [3]), unlike for a one-dimensional (ID) nonlinear os- 
cillator ,60J , as well as in contrast to the 3D cyclotron res- 
onance in a quasi-static field [56 or any wave with a par- 
allel refraction index other than unity \EV, 1551 IM] . 

The difference from Refs. ,i56t i61j i62, £3, ,641 ES] is un- 
derstood from the angle-action equations for the trans- 
verse oscillations, which now are governed by the Hamil- 
tonian 



^ = (cr - 1)5 - A V5 cos ip. 



(34) 
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FIG. 3: (color online) Particle trajectories in the momentum space (mc units) with the same po but different s: (a) po ~ 0.02, 
ao = 10~*\/2, (To = 1.008; (b), (c) po = 0, uo = 8.3, ao = 5y/2; (c) close-up of Fig. b. Blue, red, and green colors denote the 
branches 1, 2, and 3, correspondingly. All the three equilibria are center-Hke and result in stable oscillations. 



Here 3 = p^/^a is the action conjugate to the angle ip, 
and A — ao/y/a; the effective time is^,d^ = (lou/"/) dt. 
Since = (1 ~ + 2cr3)/2u [Eq. ([27|], and dr/di = 
Tp/muju, the particle motion stability in the (■0,3) space 
is equivalent to that in the 6D phase space (r,p). (We 
define -0 through x being the cylindrical phase in the 
momentum space rather than that in the coordinate 
space. In this case, Eq. (34 1 is exact; otherwise, a simi- 



lar yet approximate equation follows, in agreement with 
Refs. [SIJIS3] for ny = 1, n_L = 0, and v/c<^ 1.) 

For given a and A(ct), only one equilibrium, namely, 
a center, is allowed on the phase plane (0 mod 27r, J) 
and located aX Z = [A/2(l - ct)]^, V = (ct > 1) or 
"0 = TT (ct < 1). However, as cr = (j{ui), there exist up to 
three different phase planes, and hence the equal number 
of center points (i = 1,2,3). This situation is different 
from that, e.g., in Refs. [B2J|53], where several equilibria 
are bound to coexist on a single phase plot: as multiple 
centers are topologically impossible on a plane without 
a saddle |66], the intermediate-energy equilibrium is un- 
stable and cannot be observed there. For the topological 
constraint does not apply in our case, all the three equi- 
libria are now stable and equally realized. This results in 
unusual particle dynamics, which we discuss in Sec. |IV C| 



C. Longitudinal mass 

As the three energy states correspond to different effec- 
tive masses, a guiding center behaves differently depend- 
ing on which moff is selected; even the sign of the particle 
acceleration in response to perturbation forces can vary. 
To see this, rewrite the average motion equation (A8l as 



m,. 



dv 
dt 



(35) 



where rrin = dp/dv, or 
d 



m,, = — I -frricsv ■ _ 

ov \ 7 aw 



is the effective longitudinal mass [H71 EH] , 



d_ 
dz 



dm, 



off 



eip 



dp 
dt 



(36) 



(37) 



is the perturbation force, and p = p{z, u, t) [Eq. ( 10 )] 
A tedious yet straightforward derivation yields 



TO7' 



p3/2 



2(1 - cr)" ' 



(38) 



r2 coinciding with u'^. In the absence of the laser field 
(ao = 0), Eq. (Issl) reads m^^ = TOefT7^ > [Eq. ([l4|], as 



one would expect for a particle with meff independent of 
V [67j . Given a nonzero ao, one as well has m^^i > 0, 
because tr < 1 at the first branch, as seen from Fig. |2] 
It is also seen from Fig. |2] that > {cr()/ac)uc, where 



/l(l + 5^)1/3 [(1 + 52)1/3^(^2/2)1/3 



1/2 



(39) 



thus m||3 > 0, correspondingly. However, U2 < 
[uq/ ac)uc, yielding m||2 < for any v and s (Fig . [4j). 

Given the oscillation orbit stability (Sec. IV B[ ), a par- 
ticle residing at the second branch will exhibit unusual 
behavior in response to perturbation forces includ- 
ing gravitational and electrostatic potentials. Unlike a 
"normal" particle with a positive mass, a particle with 
m|| < will accelerate adiabatically in the direction oppo- 
site to F|| (Fig. [5]). Alternatively, should the unperturbed 
particle exhibit bounce oscillations in z (e.g., due to inho- 
mogeneity of E ot B), F^^ will shift the equilibrium point 
in the direction determined by the sign of F^Jm^^, with 
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FIG. 4: (color online) Longitudinal mass my in units m vs. 
CTo = ilo/w for a = Q, V = Q: (a) weakly relativistic case, 
flo = 10~^\/2; (b) strongly relativistic case, ao = 5\/2- The 
branches corresponding to ui, U2, and Ms are shown in blue, 
red, and green, respectively; m||i(ao = 0) = y/^ + The 
horizontal dashed lines mark zero and asymptotes at (tq — > oo: 
"1111,2 ± ao/\/2, m||3 ^ 1. The vertical dashed asymptote 
also marks the transition between the regimes with single and 
multiple branches; for conditions see the caption of Fig. |2] 
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FIG. 5: z{t) for a particle with initial u = 0, s = adia- 
batically accelerated along a magnetic field by a perturbation 
force _F|| = 10~^mao; ao = 8.3, ao = 5\/2. The sign of 
acceleration varies depending on the initial energy state 1-3 
corresponding to: (a) myi > 0, (b) m||2 < 0, (c) mya > 
(see Fig. |4|. Solid is the numerical data; dashed are analytic 
fits z{t) = F\\t^ /2m\\, with my given by Eq. (381; z and t are 
measured in units and , respectively. 



stable bouncing to persist for either sign of my (Fig. |6]). 
Merely a dissipative instability is possible for my < 
(e.g., m\\V — —rjv yields v (x e''*/l'"iil); yet it develops on 
a time scale different from that of the oscillations and, 
for weak damping, remains insignificant until large t. 

Transferring particles between the different mass 
branches also allows for a current drive effect distin- 
guishable from the traditional wave-induced methods, 
which rely on wave-induced diffusion to higher kinetic 
energies [53] • The effect is explained as follows. Sta- 



tionary fields conserve the particle quasi-energy (111 for 
a given my branch and therefore do not permit ac- 
celeration along a closed loop. However, should my 
be changed nonadiabatically along the loop, the over- 
all work performed can be nonzero; hence, even curl- 




FIG. 6: (color online) Phase plots {z,p^) showing bounce os- 
cillations along a magnetic field due to inhomogeneity of E 
or B. An additional perturbation force F\\ shifts the equilib- 
rium point in the direction determined by the sign of F\\/m\\ 
(not just F||), with stable bouncing to persist for either sign 
of my. (a,b) E{z) has a minimum at z = 0; B is uniform. 
Bounce oscillations are stable for the branches 1 (Fig. a) and 
2 (Fig. b) yet unstable for the branch 3 (particles seek high 
E\ not shown); Fy — (black) and F\\ = W~^mcuj (blue). 
(c,d) E is uniform, (c) B(z) has a maximum at 2 = 0. 
Bounce oscillations are stable for the branch 2 (shown) yet 
unstable for the branches 1 and 3 (particles seek low B; not 
shown); = (black) and F,, = 10~^mca; (blue), (d) B{z) 
has a minimum at 2 = 0. Bounce oscillations are stable for 
the branches 1 (not shown) and 3 (shown) yet unstable for 
the branch 2 (particles seek high B; not shown); same as 
in Fig. c. All figures: particles are initially at (0.1,0), s = 0; 
z and Pz are measured in units and mc, respectively. At 
the extrema: ao ~ 5\/2, and cto ~ 8.3; niyi^s > 0, m||2 < 0. 



free fields such as those due to electrostatic or pondero- 
motive potentials will be able to produce a continuous 
energy gain. Similar effects were previously discussed 
in Refs. [Tl[STJ[Mll7ni[ni[71|73J|71[75]. With the 
effective-mass formalism, these effects can now be ex- 
plained within a unified approach. 



V. CONCLUSIONS 

We showed that a classical particle oscillating in an ar- 
bitrary high-frequency or static field effectively exhibits 
a modified mass mos derived from the particle averaged 
Lagrangian [Eq. ([s])]. We obtained relativistic pondero- 
motive and diamagnetic forces, as well as magnetic drifts, 
from the rricfr dependence on the guiding center location 
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and velocity. The effective mass is not necessarily pos- 
itive and can result in backward acceleration when an 
additional perturbation force is applied. 

As an example, we explored the average motion of 
a laser-driven particle immersed in a dc magnetic field. 
Multiple energy states are realized in this case and yield 
up to three branches of rrics and the effective longitudi- 
nal mass TTT-ii for a given magnetic moment and parallel 
velocity (Fig. |4]). We showed that both mn > and 
mil < are possible then, the latter regime too allowing 
for adiabatic dynamics. From other contexts, such nega- 
tive masses are known to be capable of driving intriguing 
effects like absolute negative conductivity [751 EZ] > nega- 
tive mass instability |78l[79l[80l[8ll[82j|83l|8ail5l|86], 
and related phenomena [571 [SHI US] ■ Yet the effects that 
may flow from the variable sign of TO|| (or rricg) particu- 
larly for laser-driven particles in a magnetic field remain 
to be studied. 
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APPENDIX A: GUIDING CENTER 
LAGRANGIAN 

Consider a dynamical system, which exhibits slow 
translational motion in some guiding center variables 
(Q,P) superimposed on fast oscillations in angle-action 
variables {0,3). In the adiabatic regime, J is conserved, 
so the system action S can be put in the form 



S = 3 ■Ae+ P dQ 



Hdt, 



(Al) 



where AO is the increment of along a given trajec- 
tory, i/(Q,P;J) is the Hamiltonian, and t is the time. 
Suppose that we are only interested in guiding center tra- 
jectories, that is, those in the (Q, P) space. In this case, 
we can neglect the first term in Eq. ( |Al| ), so as to come 
up with a reduced variational principle &S = 0, where 



S 



P Q- H) dt 



(A2) 



is the new action to be varied with respect to Q and P 
only (cf., e.g., Ref. [90, §44]) 



Using S — J Ldt, where L is the Lagrangian, Eq. ( A2 ) 
can be written as |31j 



S 



L-3 e] dt 



(A3) 



By definition, the integrand here must be expressed 
in terms of the guiding center variables only. Hence 



= i/(Q,P) (a parametric dependence on J is implied 
hereafter) , so 3 ■ dt is not an exact differential, and the 
first term in Eq. (A3) is transformed as follows. By def- 
inition, L is a periodic function of 9, except that it may 
contain nonperiodic terms that are full time derivatives. 
Since omitting the latter does not affect the motion equa- 
tions, below we assume that L = (L) + Lr^, where the 
angle brackets stand for time averaging, (L^) = 0, and 
(L) is a function of (Q, P), or (Q, Q) only. On time scales 
of interest, that is. At '3> i'^^ and At 3> — {vi 
being any of the oscillation frequencies), the oscillatory 
term in Eq. (A3) vanishes; thus 



S = j (L) - 3 -u dt 



(A4) 



We can now introduce a guiding center Lagrangian 
C{Q,Q) 'AS S — J Cdt. Since the equality 



Cdt 



dt 



must hold for any time interval, one has 
£ = (L) - J • I/, 



(A5) 



(A6) 



in agreement with Refs. (91] 19^. (For a system exhibit- 
ing oscillations on multiple time scales, different Cs can 
be introduced depending on how the time averaging is 
defined.) One can also show that Eq. (A6l conforms to 
the requirement of gauge invariance: 



Rep^ 



acing L with 



L + dG/dt, where G{0, Q,t) is an arbitrary function, will 
result in C C + AC, where 



AC^u- 



do 
1)9 



Q 



dG 
dQ 



dG 



(A7) 



and G{Q,t) ^ (G(0,Q,t)). The first term is equal 
to zero due to G being periodic in 0. Thus, one has 
A£ = dG/dt, i.e., AC is, too, a full time derivative, and 
therefore does not affect the average motion equations 



d fdC \ _dC 



(AS) 



For L being a periodic function of t (rather than, or in 
addition to, 0), the above procedure would map out the 
time variable, thus yielding an analogue of the Mauper- 
tuis principle [90*, §44]. However, should t be kept (un- 
like 0) as the independent variable, the derivation of the 
guiding center Lagrangian is modified as follows. Con- 
sider the fast time t and the slow time t separately. Then 
we obtain an extended system having the action 



S^3 ■A9~HAt 



P dQ 



Hdt, 



(A9) 



where the formally introduced momentum —H conjugate 
to t is to remain constant in the adiabatic regime. The 
super-Hamiltonian H must generate the same canonical 
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equations as those of the original system; it must also 
provide that dt/dt = 1, as follows from the definition of 
t. These conditions are satisfied if one takes H — H + H, 
so the super-Lagrangian 

L = P Q + J e - Hi- H (AlO) 
equals L. Then the guiding center Lagrangian reads 

C= (L) - J -u + H, (All) 



which is equivalent to Eq. ( A6 1 , since constant H can be 
omitted. 

Results similar to those in this Appendix were ob- 
tained earlier for particle motion in a dc magnetic field 
[31 [35l IS], oscillations in nonrelativistic high-frequency 
waves [HI [S3], and laser-driven relativistic electron dy- 
namics in vacuum [16'. 1241 [^ . In the main text, we make 



use of the general form of the theorem (A6l, which con 



tains the earlier results as particular cases. This general- 
ity allows us to formulate a fundamental concept of the 
effective mass for an oscillating particle without making 
preliminary assumptions on the nature of the oscillations. 



APPENDIX B: PONDEROMOTIVE FORCES 

1. Nonrelativistic wave fields 

Let us apply the effective mass formalism to derive 
ponderomotive forces, starting with those in the nonrel- 
ativistic regime [7] [TUl [HJ [Hi [M]- Consider a particle 
oscillating in a high-frequency wave E — —Vip, where 



ip — (Po{tc, t) cos{ujt — k • r). 



(Bl) 



We will assume that kv uj; we will also assume that the 
envelope (po{r,t) varies little on the time scale u!~^ and 
has a spatial scale £ large compared to the amplitude of 
the particle oscillations {£ ^ eE/muP') and the guiding 
center displacement on the oscillation period {£ ^ v/uj) 



[7]. Then C = -? 



where 



ero 



E(f,t) 



(B2) 



is obtained using <yj(r) « — Tosc • E, with the quiver 
displacement rose = — eE/rnw'^, and the Doppler-shifted 
frequency cj' = cj — k • v. Then (£osc) = where 



4m(w — k • v)^ 



(B3) 



is known as the ponderomotive potential [7] [TOl [HJ [T5J 

IM] . Eq being the field amplitude; thus 



TTioff = m + <i>/c^ 



(B4) 



Omitting an insignificant constant, the guiding-center 
Lagrangian reads C — \ mv^ — <I>. Hence, the Hamil- 
tonian takes the well known form 



2m 



(B5) 



$ playing a role of an effective potential, as expected. 

Suppose now that, under the same conditions, an ad- 
ditional dc magnetic field B is imposed. Assuming that 

- /ii? + (Lose), where 



B is smooth, one has C = — r 
m 



M - 



2B 



(B6) 



is the new adiabatic invariant proportional to the ac- 
tion of the particle Larmor rotation at frequency f2o = 
eB/mc, Vosc is the induced oscillatory velocity propor- 
tional to E [m [MJ inS] [Ml , and (Lose) is proportional to 
E^. Suppose B « zB{z) [Eq. Eu] and take 



E = Re I 



(B7) 

where £^g"''' are smooth envelopes, and t^-'^ are unit vec- 
tors denoting polarization circular in the plane transverse 
to the dc magnetic field and that parallel to B: 



f^*^ = (A±iy)/\/2, 



.(II) 



In this case, (Lose) — 



where 



E, 



(+)|2 



E, 



(-)|2 



4mw'2 \l + no/uj' 1 - no/uj' 



E, 



(11)1 



(B8) 



(B9) 



matches the known ponderomotive potential in a dc mag- 
netic field [TOldl]. Thus, nictf = m + (fiB + ^b)/c^, and 



1 



n = —p- 

2m 



in agreement with the earlier results [55] . 



(BIO) 



2. Relativistic laser wave in vacuum 

Now consider relativistic electron motion in a vacuum 
laser field A = Ao(r,f)cos^ of arbitrary polarization, 
assuming that the vector-potential envelope Ao(r,t) has 
a scale large compared to the wavelength, ^ = uj{t—n-r/c) 
is the phase, and n = k/fc is a unit vector, say, in the z 
direction. Using Eqs. ([23|-(26l and conservation of the 



transverse canonical momentum p',^ = — (e/c)A', one has 



(a'2)^ 



(7') 



(Bll) 



where a' = eA'/mc^, and (a'^)^ = (a^)^ is an invariant. 



Using Eq. (29 1, one also has, without solving the motion 



equations, that 



(7% - 7l + (pl')e/(mc)2 ^ ^l + (a2)^, (B12) 
and £' = -mc^ (7')^. Thus, moff equals [16l [HI [H] [IH 



TTleS = m. I 



(B13) 
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which is a Lorentz invariant independent of v. Corre- 
spondingly, the guiding center momentum reads p = 
7meffV, and the well known Hamiltonian is given by 

H - ^jml^c^+p'^c^, (B14) 



the ponderomotive force resulting from the moff depen- 
dence on f and, possibly, slow dependence on t. 
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